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Abstract 

We study the most general form of a three dimensional classical integrable system with axial 
symmetry and invariant under the axis reflection. We assume that the three constants of motion 
are the Hamiltonian, H, with the standard form of a kinetic part plus a potential dependent on the 
position only, the z-component of the angular momentum, L, and a Hamiltonian-like constant, 
H, for which the kinetic part is quadratic in the momenta. We find the explicit form of these 
potentials compatible with complete integrability. The classical equations of motion, written in 
terms of two arbitrary potential functions, is separated in oblate spheroidal coordinates. The 
quantization of such systems leads to a set of two differential equations that can be presented in 
the form of spheroidal wave equations. 

1 Introduction 

This is a paper on integrability of three dimensional systems, both classical and quantum. Apart 
from being interesting by itself, this study can be extremely interesting for a field, which is rather 
far away from the present subject, the theory of the continuous media: gas, fluid and plasma. Those 
are mechanical systems with an infinite number of degrees of freedom and certainly much more 
difficult than the mechanics of a particle in three dimensions. However, in the theory of continuous 
media, we can pose a problem of the existence of soli tons: the steady solution with time-independent 
field with density p{x) and velocities v(x). It could exist in the non dissipative case only, as we 
neglect the viscosity. In this latter case, the particles which constitute the media should move in a 
self consistent potential along its trajectories, which should be closed. Although this potential is the 
result of interparticle interaction, each single particle moves in an effective potential that provides the 
closed trajectory. From the analysis of the equations of motion in continuous media (QX it follows 
that the simplest shape of solution is a toroid. Thus, we have to look for a three dimensional single 
particle integrable system including a toroidal shape. One of the goals of the present paper is to make 
a first step in this direction, as we describe all completely integrable systems with axial symmetry. 
The next step would be to select among all the completely integrable systems those sharing some 
specific properties. 
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2 The axial coordinates 



Let us consider a three dimensional classical system with the canonical coordinates (pi,Xi), i = 
1,2,3, where the commutators are defined by the usual Poisson brackets {A, B} = J^- — J|^- , 
the sum being understood over repeated indexes. The evolution of the system is described by the 
Hamiltonian 

P 2 

H = \-U(x), x = (x 1 ,x 2 ,x 3 ), p = (pi,P2,P3), (2.1) 

2m 

where U (x) is a potential term. Let us assume the existence of two additional integrals of motion for 
this system, all of them in involution. The first one will be chosen as the angular momentum along 
the direction fixed by the unit vector n = (m, "-2, TI3): 

L = n • (x x p) (2.2) 

As for the second one we require to be quadratic in momentum and commuting with both, H and L. 
Then, it is straightforward to show that, provided we include also the reflection with respect to the 
n-axis as a discrete symmetry of our system, it must take the general expression: 

H = ^- Pi g tk (x)Pk + Hx) (2-3) 
2m 

where the quadratic term has the 'metric' 

g ik (-x) = 5 ik (yL • n) 2 - x n {x ink + x k rii) + (x 2 - a 2 )n;n fc , (2.4) 
being a 2 a real constant that we take positive. The quadratic term in the momenta can be written as 

±(L 2 ± -aW), (2-5) 

where L± is the perpendicular component of L to the n-axis, and p n = p • n. Here we must point out 
that the operators (12.21) and (12.51) determine the prolate-oblate spheroidal coordinates that separate the 
Laplacian operator J91- 

The commutation of L with H and H restrict the form of their potential terms 

U(x) = C/(x 2 , (x • n) 2 ), $(x) = $(x 2 ,(x-n) 2 ) (2.6) 
while the commutation of H and H leads to the equations 

dM X )=g ik ( X )d k U( X ) (2.7) 

In order to deal with (12.71) we diagonalize the matrix g* fc (x). Its eigenvalues A(x) and eigenvectors 
-A(x) are obtained from the matrix equation 

( 5 jfc (x)-A(x)<5 ifc )^(x) = (2.8) 

We get the solutions 
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x 2 — a 2 / / x 2 — a 2 x 2 



A±(x) = ^— ±y(_^_ j +a 2 (x-n) 2 , A (x) = (x • n)* (2.9) 
The corresponding eigenvectors can be expressed in the following way: 



A+(x) = ^A_(x), 4-(x) = ^A + (x), 4 9 (x) = x2 X n)2 =frp(x) (2.10) 

where </? is the azimutal angle around the n-axis. Let us write here some useful identities of these 
eigenvalue functions for future calculations, 



(^A + ) 2 = ^(A + + a 2 ) 



(5 J A_) 2 = - x i^ r (A_ + a 2 ) (2.11) 

(^A ) 2 = -4a+A_ 
ar 

Since the eigenvectors (12.101 ) are orthogonal, it is natural to adopt as new orthogonal coordinates 
the set {A+(x), A_(x), y(x)}, which are essentially the oblate spheroidal coordinates (9). We turn 
to eq. (12.71 ) for the potentials, now expressed in the new coordinate system. Taking into account that 
due to the geometric symmetry, L (12.21 ). U and $ do not depend on tp, (12.71 ) becomes 

di\+ d + <$> + di\- = g ik {d k X+ d + U + d k X^ d^U) (2.12) 

where d± stand for Now, as g lk is diagonal in the new coordinate basis {VA_, VA + , Vy?}, this 
equation decouples in 

d + <S> = \_d + U, d^> = \ + d-U (2.13) 

or 

d+{$-\-U) = 0, d_($-A + [/) = (2.14) 

This means that 



$-A+f/ = -/(A + ), * - A_17 = -</(A_) (2.15) 

where /(A+) and <?(A_) are arbitrary functions. The minus sign in the front of ( 12.151 ) intends that the 
next equations, derived from (12.151) . can be written in a symmetric form with respect to the variables 
A + and A_. These equations give an expression for the potentials as follows: 

u.li^M, t^j^-^^y ( , 16) 

These are the most general expressions for the potentials U (x) and 3>(x) compatible with {H, H} = 
0. 
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3 Separation of variables 



3.1 The momentum as a gradient 

Now we will recall here a property for a classical system in three dimensions having three integrals 
of motion (including the Hamiltonian) hi, i = 1, 2, 3, in involution, thus being integrable: 



hi(p,x) = ei, i = 1,2,3 (3.17) 

dhidhj dhidhj . 
{hi,hj} = — —± - — —± = 0, Vz,j (3.18) 
dp k ox k ox k opk 

If we assume that the determinant of the matrix (f^l is nonvanishing, the inverse function theorem 
gives from (13.171 ) the expression of the momenta as functions of coordinates, at least locally: 

P< = /i(x,e) (3.19) 

so that 

ft i (f(x,e),x) = e< (3.20) 
Taking partial derivatives of this equation, we obtain 

dhj dfk dhj 
dpfc dx/ dx; 

d/ij d/ fc d/t m _ gfej dfr m _ dh m dhj ^ 
dpfc dx/ dp; dx; dp; dx/ dp; 

dfij d/fc d/t m | dhj dh m = Q ^ ^ 

dpfc dx/ dp; dp fc dx fc 

Note that the second identity in (13.221 ) comes after condition (13 - 1 8b - According to our hypothesis, the 
determinant of the matrix {jj^j vanishes, so that the last relation can be simplified: 



From (13.211) and (13.241 ) we have 



hence 



d fa dh m dh„, 
dxi dpi dx k 

dh m ( dh dfr \ ,„ 

' 1 " (3.25) 

(3.26) 



dpi \dxi dx k 
df k dft 



dxi dx k 

This means that the vector field, in the variable x, f (x, e) has vanishing rotational, i.e., 

rot x f (x, e) = . (3.27) 
Therefore, if f (x, e) is sufficiently regular, there exists a function F(x, e), locally defined, such that 
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flF(x,e) . e gra d x F(x,e) =f(x,e). (3.28) 

In conclusion, we have shown that if there are three integrals of motion (13.171 ) in involution, the 
momentum will be the gradient of the function F: 

dF(x,e) 

Pi = 5 (329) 

oxi 

These could be considered as a first class constraints. The function F is the characteristic function in 
the Hamilton- Jacobi approach. Note that this is a general property, valid for any dimension n. 



3.2 The separation of F 

We can apply the above results to our system by making the identification h\ = H, hi = L, h% = H. 
The next step is to show the separability of the function F in the variables {X±,ip} introduced in the 
previous section. Then, if we apply the chain rule to (13.281) and make use of (12. 10b . we have 



p = VA + 8 + F(X + , A_, ip) + VA_ <9_F(A + , A_, <p) + 2 " XX , 2 d v F(\+, \-,<p). (3.30) 

x z — (x • nj z 

Since the vector fields in (12.101) are mutually orthogonal, multiplying both sides of the above expres- 
sion by n x x, we obtain 

(nxx)-p = n- L = ^ = d !fi F(X + , \-,cp), (3.31) 

where I is the value of the integral of motion L corresponding to the angular momentum around n. 
We conclude that 

p = VA+ d + F{\ + , A_, <p) + VA_ d^F(X + , \-,<p)+ n . XX £ (3.32) 

x z — (x • nj z 

and therefore taking the square modulus in (13.321 ). we have that 

£ 2 

p 2 = (VA+) 2 (9+F) 2 + (VA_) 2 (d_F) 2 + -= — ^ = 2m (E — U) (3.33) 

x z — (x • rv) z 

where E is the value of the constant of motion H. We recall that the expression ( 13.321 ) for the mo- 
mentum vector field has been obtained in the basis that diagonalizes the metric matrix of components 
{g lk }. In this basis and from ( 13.321 ), we straightforwardly compute: 

Pi9 ik Pk = A_(VA + ) 2 (S + F) 2 + A + (VA_) 2 (9_F) 2 + ^'^^ = 2m(E - $) (3.34) 

x z — (x • n) z 

where E is the value of H. 

Next, we multiply equation (13.331) by A + , (13.341) by —1 and sum. Next we make the same proce- 
dure with A_ instead. Finally, we get the following two expressions: 
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(VA + ) 2 (A+ - X-)(d + F) 2 = 2m (\ + (E - U) - (E - $)) - ^(A + - (x • n) 2 ) 

-(VA_) 2 (A+ - A_)(cLF) 2 = 2m (\_(E - U) - (E - <&)) - 2 f (A_ - (x • n) 2 ) 

\ / x z — (x • n) z 

(3.35) 



Taking into account the definition (12.91 ) of A±, the next formulas are easily obtained: 



x 2 


-(x- 


n) 2 = 


A+ 


+ \- + a 2 


A+A_ 

+ 


= ^(A + + a 2 )(A 


A+ 


"(x 


•n) 2 = 


A + 


A+A_ _ 
■ + a* ~ 




+ a 2 ) 


A_ 


"(x 


•n) 2 = 


A_ 


A+A_ _ 


^(A + 


+ a 2 ) 


(Ah 


. - A_ 


-)(VA_ 


h) 2 


= 4A+(A+ 


+ a 2 ) 




(Ah 


.- A. 


-)(VA. 


-? 


= -4A_(A 


-+a 2 ) 





a 2 ) 



(3.36) 



Then, we carry (13.361 ) into (13.351 ) to get 

4A+(A+ + a 2 )(d+F) 2 = 2m (\ + {E -U)-(E-&) 
4A_(A_ + a 2 )(cLF) 2 = 2m (\-(E -U)-(E-$) 

We finally have 

( d+ F) 2 = J ^^i [ 2m(x + (E-U)-(E-* ) 

( ^ )2 = 4CT( 2m ( A - (£ - [/) " ( '-^ A_ + a 



A + l 2 
A+ + a 2 

A, I 2 
A_ + a 2 



A + l 2 
A + +a 2 
A_^ 2 



(3.37) 



(3.38) 



Now, let us observe that in the right hand side of (13.381 ) the combination A+ {E — U) — {E — <&) 
cannot depend on A_. At the same time, \-(E — U) — (E — $) will not depend on A+. Note that 
a similar behavior have arosen in (12- 14b - Hence, since d + F and d-F depends only on A + and A_, 
respectively, while the dependence on ip was given by (13-3 lb . we conclude that we have managed to 
obtain the separation of F in the variables A±, ip, i.e., 

F(X+,\-,<p) = A(\+)+B(\„)+£(p (3.39) 

Later, we shall find an explicit form for F in terms of a different set of variables. 

3.3 The integration of the equations of motion 

Our next objective is the integration of the following equations of motion: 

Pi = diF(x) . (3.40) 
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In other words, our goal is finding the explicit dependence of the variables A± and <p with time, i.e., 
the functions A± = \±(t) and <p = <p(t). From (13.381 ) and ( 13.391) . it is clear that 



(d+A(X + )) 2 = P(x,X+) 



(3.41) 

(cL5(A_)) 2 = P(x,A_), 
where we are using the shorthand notation 



P(x ' A *> " iMCT { 2m ( A±(£ " c,(x)> " (f " * (x)) ) " ■ (3 ' 42) 



The chain rule along p = m x and (13.401 ) give us 



1 1 

A + = <9;A + ±i = — diX + pi = — diX + diF 
m m 

A_ = di X- Xi = — diX- Pi = — diX- diF (3.43) 
m m 



ip = dup it = -5 — -7T — = — di<p diF . 

x z — (x • nj z m m 

where the dot denotes time derivative, as usual. Therefore, 



A+ = - diX+diF = - (d i X+) 2 d + A(X + ) = - — ^t- 
m m m X+ — X. 



■(X + + a 2 )d + A(X+) 



1 4A_ 



m X + — X- 



(A_ + a 2 )d_B(A_) 



a 2 l 



m (A+ + a 2 )(A_ + a 2 ) m A + -A_\A_ + a 2 A+ + 



a" 



From the first two equations in (13.44I ). 



(3.44) 



A- 



1 4(A+ + a 2 



X + d + A(X + ) m A+-A_ 

A_ _ 1 4(A_ + a 2 ) 
I X-d-B(X^) ~~m A+ - A_ 

we can get three expressions separated in A + , A_ and ip: 



(3.45) 



A- 



A_ 



X+d+A(X + ) X-d^B(X. 



4 
m 



A+(A+ + a 2 )a+^(A+) A_(A_ + a 2 )d-B{\J) 



+ 



4A+(A+ + a 2 ) 2 d + ,4(A + ) 4A_(A_ + a 2 ) 2 d_£(A 



(3.46) 
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If the potentials ( 12.161 ) are known, e.g., the functions /(A + ) and <?(A_) are chosen, then the functions 
<9_i_j4(A_|_) and d-B(X^) will be also explicitly known. Thus, integrating over time (13.461) we obtain 
the following crude expressions: 



D zd z A(z) Jo Z zd z B(z) m Cl 

A + l r x - l 

Z z(z + a*)d z A(z) + Jo dZ z(z + a?)d z B(z) ~ ° 2 ° A1) 



' Az(z + a 2 ) 2 8 z A{z) + Jo Az{z + a 2 fd z B{z) 1 + 



These relations constitute an implicit form of the integration of the equations of motion. To obtain the 
coordinates A + , A_ and if as explicit functions of time we should invert such relations, a task which 
is not often simple and that will depend on each particular choice of the functions /(A + ) and g(A_). 



3.4 The function F in terms of oblate spheroidal coordinates 

From the definition (12.91 ) A + is always positive meanwhile the values of A_ lie on the interval [—a 2 , 0] , 
depending on the scalar product x • n. Thus, we suggest the following change of coordinates: 

A+ = a 2 sinh 2 a, A_ = -a 2 sin 2 (3 . (3.48) 

Therefore we have 

x = a cosh a sin /3 cos 93, y = a cosh a sin j3 sin 93, z = a sinh a cos /? (3.49) 

where {a,/3, if} are the usual oblate spherical coordinates [Ill-O. Here, we shall see how the de- 
pendence of F in terms of a, [5 and ip gives a new insight into the above discussion. If we take time 
derivative in (13.481) . we obtain: 

A+ = (2a 2 sinh a cosh a) a ; A_ = — (2a 2 sin/3cos/3) $ . (3.50) 
From (13.481) . one readily obtains 



dF 



2 



4A + (A ++a 2)^-) = (Jl-j (3.51) 

Now, we compare the expressions given in (13.431 ) and (13.481 ) for A + and use (13.511 ) to conclude that 

• 1 1 dF 

a = rT~2 ; ■ 2 o • ( - 3 - 53 > 

m sinh a + sur oa 
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An analogous manipulation shows that 



= - 



1 



dF 



m sinh 2 a + sin 2 (3 3(3 ' 
Using the last formula in (13.441 ) and (13.481) . we find the expression for the time derivative of <p as 



(3.54) 



1 



1 



1 



(3.55) 



ma 2 sinh 2 a + sin 2 (3 \cos 2 (3 cosh a, 

Since the functions /(•) and g(-) depend respectively of A + and A_ only, they can be written as 
functions of a and (3, respectively. As /(•) and g(-) are, in principle, arbitrary, we could denote 
these functions as f(a) and g(f3) respectively. Thus, (12.151) can be written as (except for an irrelevant 
change on the sign): 



$-a 2 (sinh 2 a)U = /(a) 
Then, (13371) with (I3.51H3.52I) and ( 15351 give 



$ + a 2 (sin 2 P)U = g(p). 



(3.56) 



\da) 



8F\' 
dp) 



2m 



2m 



a 2 sinh 2 aE - E + /(a) 



sinh 2 a 
cosh 2 a 



a 2 sin 2 P E + E - g(/3) 



sin 2 /3 



cos 2 P ' 

so that, we finally get the following expressions for the derivatives of a and /?: 



(3.57) 
(3.58) 



o 



ft 



1 



a 2 m sinh 2 a + sin 2 P 
1 1 



2m 



aV 2 



a 2 sinh 2 a£-B + /(a) 



sinh 2 a 
cosh 2 a 



1/2 



a 2 m sinh 2 a + sin 2 P 
1 1 



a 2 m sinh 2 a + sin 2 /? 



2m 



a 2 sin 2 PE + E~g(p) 



sin 2 g 
cos 2 /9 



.1/2 



a^m sinh a + sin" 1 /? 



1/2 



where obviously, 



(3.59) 



(3.60) 



A, 



2m 



2m 



a 2 sinh 2 a E — E + /(a 
a 2 sin 2 ^£ + 



„ 9 sinh a 
cosh a 

sin 2 /? 
cos 2 P 



(3.61) 
(3.62) 
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Clearly, after (15391 and (1X601 . we have 



a 







A V2 A l/2 







cia 



o A 



1/2 



d/3 



o A 



1/2 



c 2 , 



(3.63) 



where C2 is a constant with respect to time (obviously, it depends on a and j3). Another constant of 
motion can be obtained as follows: First, we write (13.551 ) as 



1 



1/9 ma? sinh 2 a + sin 2 (3 I cos 2 (3 
Then, using (15391 and (13351 in (I3T64T) . we have 



cos 2 /3 + sin 2 /3 cosh 2 a — sinh 2 a 



cosh a 



(3.64) 



From (13.651) . we obtain 

d_ 

It 

which shows that 



P sin 2 



+ 



sinh 2 a 

A*/ 2 cos 2 (3 ' A V 2 cosh 2 a 



<ia sinh 2 a f d(3 sin 2 /3 



a 2 w a^ 2 cosh 5 



0, 



(3.65) 



(3.66) 



¥ ~ — 



da sinh 2 a f d(3 sin 2 /3 

7V2" COS 2 /? 

/3 



+ 



c 3 



A^ 2 cosh 2 a 
is a constant of motion. 

Furthermore, from ( 13.591 ) and ( 13.601 ) we can obtain a third result as we show in the sequel 



(3.67) 



a 



(3 . 2 1 sinh 2 a 1 sin 2 /? 1 



—- sinh 2 a H m sin 2 /3 = . : 2""' . 2 3 + "T ■ 1 2 , ■ 2 a = ~ ' ( 3 " 68) 

A„' Ai' a m sinh a + sin p arm sinh a + sin p m 

which obviously yields after integration 

f da 2 f d{3 . o n 1 

/ — -7^ smh a + / — -7^ sm p = — i . (3.69) 

J a« 7 •/ Aj a m 

The function F can now be written in terms of the variables (a, (3, ip). Since A + and A_ are functions 
of a and (3 alone, respectively, formula (13.391 ) can be written as 

F{a,p,<p)=A{a) + B{p)+l<p (3.70) 

where 

r dF f dF 

A(a)= —da, B((3)= —d@ (3.71) 



da J dp 
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This gives the final expression for F(a, (3, cp) as 



F(a, P,<p) = J A l J 2 da + J A X J 2 d(3 + £<p. (3.72) 
Time invariants C2 and C3 can be written in terms of certain partial derivatives of F(a>, f3, tp) as we 



can easily show. In fact, using the expressions for A]/ 2 and A l J 2 in (13.611) and (13.621 ). we have that 



9F 2 „ dF „ dF 

= - m 2 C 2 , -^7 = ^3, 7^=^ ( 3 - 73 ) 

as it can be easily checked. Then, the function F can be written as 

F = Et — mEC 2 + £C 3 , (3.74) 
where C2 and C3 are dependent on a, (3 and ip, but they are time independent. 



4 Quantum systems 

From the point of view of quantum mechanics, the Hamiltonian H as well as the integrals of motion 
L and H are hermitian operators obtained simply by replacing pj, —>■ —i dp., k = 1,2, 3, in ( 12. IK (12.21 ) 
and (I2.3I ). respectively (we have taken ft = 1 along this section). Formal hermiticity follows from 
the fact that the operators H, L and H are symmetric in the usual cartesian coordinates. The kinetic 
parts of H and H are nonsingular quadratic expressions on positions and momenta and they are self 
adjoint ( ifTTI ). In addition, we assume that the potentials U(x) and 3>(x) satisfy sufficient conditions 
so that both H and H be self adjoint (as for example that the conditions in the Kato Rellich theorem 
be satisfied lfl2l ). 

The conditions (12.41 ) on the metric g(x), and (12.6I )- (12.7I ) on the potential terms U (x), $(x) guar- 
antee the commutation relations for these operators: 

[H,L) = [H,L] = [H,H]=0. (4.75) 

We will look for the simultaneous eigenf unctions ip(x) of the three operators 

(H - E)ip = (H- E)ip = (L - £)tp = (4.76) 

or equivalently 

-AV>(x) = 2m{E - U(x))tp(x) 

-A<0(x) = 2m(E - $(x))^(x) , (4.77) 
-idtpipfa) = £-0( x ) 

where 

A = d k d k , A = d j{ f k {x)d k . (4.78) 
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Now, we will express these differential operators in terms of the coordinates A + , A_ , ip in order to 
rewrite the eigen-equations in the form 



A+ - A_ 



4A+(A+ + a 2 )V>++ + 2(a 2 + 3A+)^ H 
4A_(A_ + a 2 )^__ + 2(a 2 + 3A_)V>_ - 



X + + a^ w 



A_ + a 2 ^ w 



: 2m(E - U)ip 



A_ 



A+ - A, 



4A+(A+ + a 2 )V>++ + 2(a 2 + 3A + )^+ + 



A. 



+ 



A+-A_ 



X + + a^ w 
4A+(A+ + a 2 )V>— + 2(a 2 + 3A_)V>- + 



A_ + a 2 ' /W 
Taking into account that 

and also (12.141) . equations ( 14.791 ) can be written in a separated form, 



2m(E - 

(4.79) 
(4.80) 



4A + (A + + a 2 )V>++ + 2(a 2 + 3A+)V>+ - 



A+ + a 2 

4A_(A_ + a 2 )V— + 2(a 2 + 3A_)V- - f A ~ 9 



2m 
2m 



(E - U)X+ - (E — $)] V 



(£ - £7)A_ -(E- $)] V 



where A+ and A_ appear in the first and second equation in (14.811 ) respectively. Therefore 
look for a factorized solution for the eigenf unctions as follows: 



(4.81) 
we can 



VKx)=V + (A + )V-(A_)e^ (4.82) 

We can go back to the change of coordinates given by (13.48I ). This change of coordinates is also 
suggested by the formulas below, as we shall see. Now, it is time for choosing explicit forms for the 
functions f(x) and g(x) in (12.151) . For f(x), we shall choose the function that vanish identically. For 
g(x), we choose <?(A_) := — Q(\- + a 2 ), where Q is a constant. Then, the following expressions 
arise: 



/(A+) = $ - X+U = 0, g(X-) = $ - X-U = -Q(A_ + a 2 



C/(x) = -Q 



Or, in cartesian coordinates, 



A_ + q 2 
A+ - A_ ' 



<f(x) = -Q 



x 2 + a 2 



A+(A_+a 2 ) 
A+ - A_ 



2 V[(x 2 + a 2 ) 2 -4a 2 (x 2 -(x-n) 2 )] 



2MV2 



1 



(4.83) 



(4.84) 



Now, we carry ( 13.481 ) and ( 14.831 ) into ( 14.811 ) to get the following set of two differential equations in 
the variables a and (3: 
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d 2 sinh a d sinh 2 a 
da 2 cosh a da cosh 2 a 



+ 2m 



E a 2 sinh 2 a — E 



ip + (a) = 0, (4.85) 



d 2 sin/3 cZ sin 2 /? 
<i/3 2 cos (3 dp cos 2 /3 



' 2 + 2m 



^a^sin^/^ + .E + Qa^cos 2 /? ^~(^) = 0. (4.86) 



This is a special type of equations already studied in the literature that we briefly analyze in the next 
section. 
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Figure 1: Iso-potential lines in the XZ-plane for three values of the parameters. 



5 Study of the equations and their solutions. 



First of all, it seems convenient to simplify equations (14. 85114. 86b - In order to fulfil this goal, let us 
choose the following new coordinates: 
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t = sinha, u = sin/3, - — = \J t 2 + 1 — , — — = \J\ — u 2 — . (5.87) 

da at dp du 

Let us introduce the following parameters: 

£ = 2ma 2 E, q = 2ma 2 Q, £ = £ 2 + 2mE , (5.88) 
then, eq. (14.85H4.86l ) become respectively 



/2 ,d 2 d £ 2 



+ £t 2 -£ 



= 0. 



(5.89) 



(1 



^ (in 2 cfu 



1 



£u 2 + o(l - u 2 ) + £ 



=0. 



(5.90) 



We observe that in the case of £ = q = and only in this case, these equations can be reduced 
to equations of hypergeometric type, which can be solved in terms of hypergeometric functions. 
However, this is not the most general case, let us consider the following differential equation: 



(1 " z 2 ) 



d 2 w 
H 2 



dw [ , 
~dz + I 



1 



+ 7 2 (l-z 2 Uu> = 0, 



(5.91) 



where A, u and j 2 are real parameters (7 2 may be positive or negative depending on the use of prolate 
or oblate coordinates respectively). This is called the spheroidal wave equation ifTOl . It is very simple 
to show that both equations (15.891 ) and (15.901) are versions of the spheroidal wave equation. In fact, 
( 15.901 ) can be written as 



(1 - u 2 )^^ - 2u ^-M + j G - + ,'(1 - u 2 ) } r in) = 0, (5.92) 

with G := £ + £ and q' = q — £. If we apply the change of variables given by t = ia, equation 
( 15.891 ) becomes 

^(Ptb+ia) dib+(a) { t 2 , 9 , 1 , , , 

where again G = £ + £ and we have kept the notation t/j + (a) = t/j + (ia) = i/) + (t). 

Solutions of the spheroidal wave equation ( 15.911 ) and therefore of ( 15.921 ) and ( 15.93b have been 
studied in |[T0l . The origin z = is a regular point of the equation and therefore, we can find two 
linearly independent functions in terms of powers series on the variable z. These series converges 
on the open circle centered at the origin and radius equal to one, since z = ±1 are singular points 
for the equation. On this open circle, one can find one even and one odd solution of (15.911 ) of the 
form J2^=o a n z 2n and J2^=o a n z 2n+l respectively, which are linearly independent. These series 
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do not converge at the singular points ±1. As they do converge on the open interval (—1,1), the 
wave function ip~ solution of equation (14.861) is periodic on the real axis with singularities at the 
points (2n + l)7r/2. There exists another type of linearly independent even and odd solutions on a 
neighborhood of the origin that may converge at the singular points ±1, provided that a relation is 
satisfied between the coefficients A, fi and 7 in (15.911 ) (and its corresponding translation in terms of 
the coefficients in ( 15.92b and ( 15.931 )) flOl . In any case, the solutions tp~(f3) of (14.86! ) on the real axis 
are periodic and therefore, not square integrable. 

With respect to equation (15.91I ). z = ±1 are regular singular points with indices equal to ±^yu = 
±i£. Being I an integer, the two linearly independent solutions on a neighborhood of z = 1 are u\{z) 
and u 2 (z) with u x (z) = (z - l) £ / 2 Fi(z - 1) and u 2 (z) = u x (z) log(z -l) + (z- l)~ l / 2 F 2 {z - 1), 
where F\(z—1) and F 2 (z — 1) are power series on z — 1. These power series have radii of convergence 
equal to 2. On a neighborhood of z = — 1, similar solutions can be found. Power series never truncate. 

There is another singular point at z = 00. This singular point is irregular. Solutions on a neigh- 
borhood of the infinite have the form z v Y^ = _ OQ a nZ 2n , where v is a complex number depending 
on the equation parameters. In order to simplify the recurrence relations for the coefficients a n , it is 
customary to choose this solution as (z 2 — l)^/ 2 z u ~^ Ylri=-aa a nZ 2n ■ The condition that the Laurent 
series converges in 1 < \z\ < 00 gives a relation between A, fi, 7 and v [10]. These solutions are of 
the form 



00 

{z 2 -iyl 2 z» <n4+2nM, j = 1,2,3,4, (5.94) 

n=— oo 

where 



with zP(z) 



J u (z), Zi 2 \z) = Y v {z), Zf\z) 



Hl 1] (z) and zi A \z) 



(5.95) 
h1 2 \z), being 



(i) 

J v (z),Y v (z) and H v (z) the Bessel functions of first, second and third class respectively. Any two 
of the set of solutions ( 15.941 ) are linearly independent provided that v be not a half odd integer. 

Solutions of (15.891 ) and ( 15.901 ) can be obtained without resorting to the standard study of the 
spheroidal wave function. For instance, if we use the change of variables given by z := t 2 in (15.891 ), 
this equation is transformed into 



4^ + l,^ + (6z + 2)A + _^_ 



+ £z-£ 



0. 



(5.96) 



Now, the singular regular points lie at z = 0,-1 and it is not difficult to obtain solutions in form 
of power series on a neighborhood of these points. For example, for z = the characteristic expo- 
nents are and 1/2 giving respective linearly independent solutions of (15.961 ) of the form ipo(z) = 

J2^=o a n z " an d tpi/ 2 ( z ) = J2^=obn zn+ ^ ■ Recurrence relations for the coefficients depend on four 
coefficients, except the first and second relations which depend on the two and three first coefficients 
respectively (which is compatible with the fact that ao and 60 should be the only independent coeffi- 
cients). On a neighborhood of the singular point z = — 1, two linearly independent solutions can be 
found of the form xjj x (z) = E^=o°n(^ + l) n+£/2 and ^ 2 {z) = ij) X {z) \og{z + 1) + En=o b n(z + 
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l) n ^/ 2 . These series make sense provided that compatibility relations exists between the parameters 
t, £ and £ in complete agreement with the general study of the solutions of spheroidal wave functions 

iniia. 

6 Conclusions and remarks. 

We have studied the conditions of integrability of a classical or quantum system having a symmetry 
axis. As in a three dimensional integrable system, we have found three independent observables such 
that their respective Poisson brackets are zero, in the classical case, or commute in the quantum case. 
The chosen symmetry forces one of the observables to be the component of the angular momentum 
in the direction of the symmetry axis. The other two can be written in Hamiltonian form as a sum of 
a kinetic term plus a potential. 

In the classical case, we have obtained the most general form of the potentials corresponding to 
both Hamiltonians in terms of oblate spheroidal coordinates, that depends on two arbitrary functions 
depending on one coordinate only. We have written the equations of motion in terms of this coor- 
dinates and show that the Hamilton- Jacobi characteristic function can be written as a sum of three 
functions each one depending on one coordinate only. Then, we have obtained the explicit form for 
these three functions. 

The quantum case is obtained by direct canonical quantization of the classical case. The condition 
of integrability yields to two Schrodiger type equations in which with separate variables. Then, a 
reasonable choice on the functions that determine the potentials yields to new equations that are 
shown to be of the spheroidal type. We finish the discussion with some comments on the solutions of 
this kind of equations. 
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